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Abstract 



' The X(j)^ model in a finite volume is studied within a non-gaussian Hartree- 

■ Fock approximation (tdHF) both at equilibrium and out of equilibrium, with 

particular attention to the structure of the ground state and of certain dy- 
I namical features in the broken symmetry phase. The mean-field coupled 

I time-dependent Schroedinger equations for the modes of the scalar field are 

' derived and the suitable procedure to renormalize them is outlined. A fur- 

Q_i' ther controlled gaussian approximation of our tdHF approach is used in order 

^ . to study the dynamical evolution of the system from non-equilibrium initial 

D i conditions characterized by a uniform condensate. We find that, during the 

' slow rolling down, the long-wavelength quantum fiuctuations do not grow to a 

macroscopic size but do scale with the linear size of the system, in accordance 
^ ■ with similar results valid for the large N approximation of the 0{N) model. 

^ . This behavior undermines in a precise way the gaussian approximation within 

our tdHF approach, which therefore appears as a viable mean to correct an 
unlikely feature of the standard HF factorization scheme, such as the so-called 
"stopping at the spinodal line" of the quantum fiuctuations. We also study 
the dynamics of the system in infinite volume with particular attention to the 
asymptotic evolution in the broken symmetry phase. We are able to show 
that the fixed points of the evolution cover at most the classically metastable 
part of the static effective potential. 
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I. INTRODUCTION 



A great effort has been devoted in the last few years in order to develop a deeper qualita- 
tive and quantitative understanding of systems described by interacting quantum fields out 
of equilibrium. There is a class of physical problems that requires the consistent treatment 
of time dependent mean-fields in interaction with their own quantum or thermal fluctua- 
tions. We may mention, among others, the problem of reheating of the universe after the 
inflationary era of exponential growth and cooling, and the time evolution of the scalar order 
parameter through the chiral phase transition, soon to be probed in the forthcoming heavy- 
ion experiments at CERN-SPS, BNL-RHIC and CERN-LHC. In these situations, a detailed 
description of the time-dependent dynamics is necessary to calculate the non-equilibrium 
properties of the system. Indeed, the development of practical general techniques and the 
advent of faster and cheaper computers have made possible the discovery of novel and un- 
expected phenomena, ranging from dissipative processes via particle production to novel 
aspects of symmetry breaking 

From the technical point of view, it should be pointed out, first of all, that a perturbative 
treatment of this dynamical problem is meaningful only when the early time evolution is 
considered. The presence of parametric resonant bands or spinodal instabilities (in the case, 
respectively, of unbroken or spontaneously broken symmetries) rapidly turns the dynamics 
completely non-linear and non-perturbative. Thus, the asymptotic evolution at late time 
can be consistently studied only if approximate non-perturbative methods are applied to 
the problem 

Quite recently one of these schemes, namely the large expansion at leading order 
(LN) has been used in order to clarify some dynamical aspects of the 0^ theory in 

3 spatial dimensions, reaching the conclusion that the non-perturbative and non-linear 
evolution of the system might eventually produce the onset of a form of non-equilibrium 
Bose-Einstein condensation (EEC) of the long-wavelength Goldstone bosons usually present 
in the broken symmetry phase [0,^,0. Another very interesting result in [|^ concerns the 
dynamical Maxwell construction, which reproduces the fiat region of the effective potential 
in case of broken symmetry as asymptotic fixed points of the background evolution. 

In a companion work we have addressed the question of whether a standard BEG 
could actually take place as time goes on, by putting the system in a finite volume (a 
periodic box of size L) and carefully studying the volume dependence of out-of-equilibrium 
features in the broken symmetry phase. We summarize here the main result contained in 
PJ. The numerical solution shows the presence of a time scale tl, proportional to the linear 
size L of the system, at which finite volume effects start to manifest, with the remarkable 
consequence that the zero-mode quantum fluctuations cannot grow macroscopically large if 
they start with microscopic initial conditions. In fact, the size of low-lying widths at time 
Tl is of order L, to be compared to order L^/^ for the case of standard BEG. In other words 
we confirmed that the linear growth of the zero mode width, as found also by the authors of 
0,^]^, really signals the onset of a novel form of dynamical BEG, quite different from the 
standard one described by equilibrium finite-temperature field theory. This interpretation 
is reinforced by the characteristics of the long- wavelength fluctuations' spectrum. 

Since after all the large A^ approximation is equivalent to a Gaussian ansatz for the 
time-dependent density matrix of the system p|J^, one might still envisage a scenario in 
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which, while gaussian fluctuations would stay microscopic, non-gaussian fluctuations would 
grow in time to a macroscopic size, leading to an occupation number for the zero mode 
proportional to the volume Lr" of the system. Therefore, in order to go beyond the gaussian 
approximation, we will consider in this work a time-dependent HF approach capable in 
principle of describing the dynamics of some non-gaussian fluctuations of a single scalar 
fleld with 0^ interaction. 

Before going into the details of the analysis, let us briefly summarize the main limi- 
tations and the most remarkable results of the study of a scalar fleld out of equilibrium 
within the gaussian HF scheme |]T],|TD|-|T^] . First of all, this scheme has the advantage of 
going beyond perturbation theory, in the sense that the (numerical) solution of the evolu- 
tion equations will contain arbitrary powers of the coupling constant, corresponding to a 
non-trivial resummation of the perturbative series. For this reason, the method is able to 
take into account the quantum back-reaction on the fluctuations themselves, which shuts 
off their early exponential growth. This is achieved by the standard HF factorization of 
the quartic interaction, yielding a time dependent self-consistently determined mass term, 
which stabilizes the modes perturbatively unstable. The detailed numerical solution of the 
resulting dynamical equations clearly shows the dissipation associated with particle produc- 
tion, as a result of either parametric ampliflcation in case of unbroken symmetry or spinodal 
instabilities in case of broken symmetry, as well as the shut off mechanism outlined above. 

However, the standard HF method is really not controllable in the case of a single scalar 
fleld, while it becomes exact only in the N oo limit. Moreover, previous approaches 
to the dynamics in this approximation scheme had the unlikely feature of maintaining a 
weak (logarithmic) cut-off dependence on the renormalized equations of motion of the order 
parameter and the mode functions 

In this article we consider the case of a single scalar fleld (i.e. = 1). With the aim of 
studying the dynamics of the model with the inclusion of some non-gaussian contributions, 
we introduce an improved time-dependent Hartree-Fock approach. Even if it is still based on 
a factorized trial wavefunction(al), it has the merit to keep the quartic interaction diagonal 
in momentum space, explicitly in the hamiltonians governing the evolution of each mode of 
the fleld. In this framework, issues like the static spontaneous symmetry breaking can be 
better understood, and the further gaussian approximation needed to study the dynamics 
can be better controlled. In particular, questions like out-of-equilibrium "quantum phase 
ordering" and "dynamical Bose-Einstein condensation" can be properly posed and answered 
within a veriflable approximation. 

We also perform a detailed study of the asymptotic dynamics in inflnite volume, with 
the aim of clarifying the issue of Maxwell construction in this approximation scheme. In 
fact, in the 0{N) $^ model at leading order, the asymptotic dynamical evolution of the 
mean fleld completely covers the spinodal region of the classical potential, which coincides 
with the flatness region of the effective potential. This is what is called dynamical Maxwell 
construction 0. When we use the HF approximation for the case of iV = 1, we flnd that 
the spinodal region and the flatness region are different and the question arise of whether a 
full or partial dynamical Maxwell construction still takes place. 

In section ^ we set up the model in flnite volume, deflning all the relevant notations and 
the quantum representation we will be using to study the evolution of the system. 

We introduce in section |T| our improved time-dependent Hartree-Fock (tdHF) ap- 
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proximation, which generahzes the standard gaussian self-consistent approach [0 to 
non-gaussian wave-functionals; we then derive the mean-field coupled time-dependent 
Schroedinger equations for the modes of the scalar field, under the assumption of a uniform 
condensate, see eqs ( p.5|) , (|3.6| ) and ( p.7|) . A significant difference with respect to previous 
tdHF approaches |^ concerns the renormalization of ultraviolet divergences. In fact, by 
means of a single substitution of the bare coupling constant Ab with the renormalized one 
A in the Hartree-Fock hamiltonian, we obtain cut-off independent equations (apart from 
corrections in inverse powers, which are there due to the Landau pole). The substitution is 
introduced by hand, but is justified by simple diagrammatic considerations. 

One advantage of not restricting a priori the self-consistent HP approximation to gaussian 
wave-functionals, is in the possibility of a better description of the vacuum structure in case 
of broken symmetry. In fact we can show quite explicitly that, in any finite volume, in 
the ground state the zero-mode of field is concentrated around the two vacua of the 
broken symmetry, driving the probability distribution for any sufficiently wide smearing of 
the field into a two peaks shape. This is indeed what one would intuitively expect in case 
of symmetry breaking. On the other hand none of this appears in a dynamical evolution 
that starts from a distribution localized around a single value of the field in the spinodal 
region, confirming what already seen in the large approach More precisely, within a 
further controlled gaussian approximation of our tdHF approach, one observes that initially 
microscopic quantum fiuctuations never becomes macroscopic, suggesting that also non- 
gaussian fiuctuations cannot reach macroscopic sizes. As a simple confirmation of this fact, 
consider the completely symmetric initial conditions (0) = (0) = for the background: in 
this case we find that the dynamical equations for initially gaussian field fiuctuations are 
identical to those of large A^ (apart for a rescaling of the coupling constant by a factor of 
three; cfr. ref. [§]), so that we observe the same asymptotic vanishing of the effective mass. 
However, this time no interpretation in terms of Goldstone theorem is possible, since the 
broken symmetry is discrete; rather, if the width of the zero-mode were allowed to evolve 
into a macroscopic size, then the effective mass would tend to a positive value, since the 
mass in case of discrete symmetry breaking is indeed larger than zero. 

Anyway, also in the gaussian HP approach, we do find a whole class of cases which exhibit 
the time scale r^. At that time, finite volume effects start to manifest and the size of the 
low-lying widths is of order L. We then discuss why this undermines the self-consistency 
of the gaussian approximation, imposing the need of further study, both analytical and 
numerical. 

In section |^ we study the asymptotic evolution in the broken symmetry phase, in 
infinite volume, when the expectation value starts within the region between the two minima 
of the potential. We are able to show by precise numerical simulations, that the fixed 
points of the background evolution do not cover the static fiat region completely. On the 
contrary, the spinodal region seems to be absolutely forbidden for the late time values of 
the mean field. Thus, as far as the asymptotic evolution is concerned, our numerical results 
lead to the following conclusions. We can distinguish the points lying between the two 
minima in a fashion reminiscent of the static classification: first, the values satisfying the 
property w/a/S < 0oo < v are metastable points, in the sense that they are fixed points 
of the background evolution, no matter which initial condition comprised in the interval 
(— we choose for the expectation value 0; secondly, the points included in the interval 
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< 0OO < are unstable points, because if the mean field starts from one of them, 

after an early slow rolling down, it starts to oscillate with decreasing amplitude around 
a point inside the classical metastable interval. Obviously, = w is the point of stable 
equilibrium, and = is a point of unstable equilibrium. Actually, it should be noted that 
our data do not allow a precise determination of the border between the dynamical unstable 
and metastable regions; thus, the number we give here should be looked at as an educated 
guess inspired by the analogous static classification and based on considerations about the 
solutions of the gap equation [see eq. ( |3.36| )] 



Finally, in section we give a brief summary of the results presented in this article 
and we outline some interesting open questions that need more work before being answered 
properly. 

II. CUTOFF FIELD THEORY 

Let us consider the scalar field operator and its canonically conjugated momentum tt 
in a D— dimensional periodic box of size L and write their Fourier expansion as customary 

0(x) = 0, e^^-^ , 0l = 

k 



D 



with the wavevectors k naturally quantized: k = {2TT/L)n, n G Z 

The canonical commutation rules are [0^ , T^-k'] = "i^kkK as usual. The introduction of a 
finite volume should be regarded as a regularization of the infrared properties of the model, 
which allows to "count" the different field modes and is needed especially in the case of 
broken symmetry. 

To keep control also on the ultraviolet behavior and manage to handle the renormalization 
procedure properly, we restrict the sums over wavevectors to the points lying within the 

dimensional sphere of radius A, that is /c^ < A^, with Af = AL/27r some large integer. 
Till both the cut-offs remain finite, we have reduced the original field-theoretical problem to 
a quantum-mechanical framework with finitely many (of order Af^^^) degrees of freedom. 

The 0^ Hamiltonian is 



7r2 + (a0)^ + m2 02 + Ab0'] 



A 



where and Ab are the bare parameters and depend on the UV cutoff A in such a way to 
guarantee a finite limit A ^ oo for all observable quantities. It should be noted here that, 
being the theory trivial (as is manifest in the resummed one-loop approximation due 
to the Landau pole) the ultraviolet cut-off should be kept finite and much smaller than the 
renormalon singularity. In this case, we must regard the 0^^ model as an effective low-energy 
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theory (here low-energy means practically all energies below Planck's scale, due to the large 
value of the Landau pole for renormalized coupling constants of order one or less). 
We shall work in the wavefunction representation where {(fl \E') = "^{(f) and 

d 

O(po 

while for /c > (in lexicographic sense) 

1 1 f d d \ 

{(p±k'^){^) = -j=iifik±i V-k) "^i^) , (vr±fc^)(^) = I ± ^(V') 

Notice that by construction the variables ipk are all real. 

In practice, the problem of studying the dynamics of the 0^ field out of equilibrium 
consists now in trying to solve the time- dependent Schroedinger equation given an initial 
wavefunction \l/(v?, t = 0) that describes a state of the field far away from the vacuum. This 
approach could be very well generalized in a straightforward way to mixtures described by 



density matrices, as done, for instance, in [10,15,16|. Here we shall restrict to pure states. 



for sake of simplicity and because all relevant aspects of the problem are already present in 
this case. 

We shall consider here the time-dependent Hartree-Fock (tdHF) approach (an improved 
version with respect to what is presented, for instance, in [0), being the large N expansion 
to leading order treated in another work [§. In fact these two methods are very closely 
related (see, for instance in [1^). However, before passing to any approximation, we would 
like to stress that the following rigorous result can be immediately established in this model 
with both UV and IR cutoffs. 



A. A rigorous result: the effective potential is convex 

This is a well known fact in statistical mechanics, being directly related to stability 
requirements. It would therefore hold also for the field theory in the Euclidean functional 
formulation. In our quantum-mechanical context we may proceed as follow. Suppose the 
field (j) is coupled to a uniform external source J. Then the ground state energy Eq{J) is 
a concave function of J, as can be inferred from the negativity of the second order term 
in AJ of perturbation around any chosen value of J. Moreover, Eq{J) is analytic in a 
finite neighborhood of J = 0, since J(f) is a perturbation "small" compared to the quadratic 
and quartic terms of the Hamiltonian. As a consequence, this effective potential Kff(0) = 
Eq{J) — J(j), (j) = E'q{J) = {(j))Q, that is the Legendre transform of Eq{J), is a, convex analytic 
function in a finite neighborhood of = 0. In the infrared limit L oo, Eq{J) might 
develop a singularity in J = and Kfr(0) niight flatten around = 0. Of course this 
possibility would apply in case of spontaneous symmetry breaking, that is for a double-well 
classical potential. This is a subtle and important point that will play a crucial role later on, 
even if the effective potential is relevant for the static properties of the model rather than 
the dynamical evolution out of equilibrium that interests us here. In fact such evolution is 
governed by the CTP effective action ||18|,19| and one might expect that, although non-local 



in time, it asymptotically reduces to a multiple of the effective potential for trajectories 
of (f){t) with a fixed point at infinite time. In such case there should exist a one-to-one 
correspondence between fixed points and minima of the effective potential. 
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III. TIME-DEPENDENT HARTREE FOCK 



In order to follow the time evolution of the non-gaussian quantum fluctuations we con- 
sider in this section a time-dependent HF approximation capable in principle of describing 
the dynamics of non-gaussian fluctuations of a single scalar field with 0^ interaction. 

We examine in this work only states in which the scalar field has a uniform, albeit 
possibly time-dependent expectation value. In a tdHF approach we may then start from a 
wavefuction of the factorized form (which would be exact for free fields) 



k>0 



(3.1) 



The dependence of ipk on its two arguments cannot be assumed to factorize in general 
since space translations act as 5*0(2) rotations on (pk and (p-k (hence in case of translation 
invariance ^pk depends only on fk~^dk)- The approximation consists in assuming this form 
as valid at all times and imposing the stationarity condition on the action 



6 Jdt{zdt-H)=0, (.) = (^(t)|.|^(t)) 



(3.2) 



with respect to variations of the functions ipk- To enforce a uniform expectation value of 
we should add a Lagrange multiplier term linear in the single modes expectations {(fk) for 
k ^ 0. The multiplier is then fixed at the end to obtain (ifk) = for all k ^ 0. Actually one 
may verify that this is equivalent to the simpler approach in which (ipk) is set to vanish for all 
k ^ before any variation. Then the only non trivial expectation value in the Hamiltonian, 
namely that of the quartic term, assumes the form 



2\2 



2LD 



T.kd+dkf)-^{(.d)+{d-k)y 

k>0 '- 



(3.3) 



Notice that the terms in the first row would cancel completely out for gaussian wavefunctions 
ipk with zero mean value. The last term, where the sum extends to all wavevectors k, 
corresponds instead to the standard mean field replacement (0^) — > 3(0^)^. The total 
energy of our trial state now reads 



E 



1 



Ab 



(3.4) 



and from the variational principle ( |3.2| ) we obtain a set of simple Schroedinger equations 

idtipk = Hk^k (3.5) 



Ho 
Hk 



15' 1 2 2 

■2 9^ + 2"°^° + 4L^ 



Ab 4 



M 



+ 



3Ab 



{d + d-k) 



(3.6) 
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which are coupled in a mean-field way only through 



1 

Id 



(3.7) 



<j2<a2 
q^k.-k 



and define the HF time evolution for the theory. By construction this evolution conserves 
the total energy E of eq. (p.4|). 

It should be stressed that in this particular tdHF approximation, beside the mean-field 
back-reaction term of all other modes on ujI, we keep also the contribution of the diagonal 
scattering through the diagonal quartic terms. In fact this is why has no contribution 
from the fc— mode itself: in a gaussian approximation for the trial wavefunctions ipk the 
Hamiltonians if^ would turn out to be harmonic, the quartic terms being absent in favor of 
a complete back-reaction 



1 



k 



(3.8) 



Of course the quartic self-interaction of the modes as well as the difference between E 
and Efc are suppressed by a volume effect and could be neglected in the infrared limit, 
provided all wavefunctions stays concentrated on mode amplitudes of order smaller 
than L^/^. This is the typical situation when all modes remain microscopic and the volume 
in the denominators is compensated only through the summation over a number of modes 
proportional to the volume itself, so that in the limit L — oo sums are replaced by integrals 



Indeed we shall apply this picture to all modes with k ^ 0, while we do expect exceptions 
for the zero-mode wavefunction ipQ. 

The treatment of ultraviolet divergences requires particular care, since the HF approxi- 
mation typically messes things up (see, for instance, |^). Following the same login of the 
large N approximation ||l]J^J^, we could take as renormalization condition the requirement 
that the frequencies ul are independent of A, assuming that ml and Ab are functions of A 
itself and of renormalized A— independent parameters and A such that 



k^ + + 3A [Sfc] 



finite 



(3.9) 



where by [.jfinite "we mean the (scheme-dependent) finite part of some possibly ultraviolet 
divergent quantity. Unfortunately this would not be enough to make the spectrum of energy 
differences cutoff-independent, because of the bare coupling constant Ab in front of the 
quartic terms in H^. and the difference between E and E^ [such problem does not exist 
in large because that is a purely gaussian approximation]. Again this would not be a 
problem whenever these terms become negligible as L — oo. At any rate, to be ready to 
handle the cases when this is not actually true and to define an ultraviolet-finite model also 
at finite volume, we shall by hand modify eq. (|3.3|) as follows: 



Ab / d^x (0(a:)^) =XL-D \ {^D 



-D 



fc>0 



+ 3Ab L^T? 



(3.10) 
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We keep the bare coupling constant in front of the term containing because that part 



of the hamiltonian is properly renormalized by means of the usual cactus resummation ||2T 



which corresponds to the standard HF approximation. On the other hand, within the same 
approximation, it is not possible to renormalize the part in curly brackets of the equation 
above, because of the factorized form ( |3.1| ) that we have assumed for the wavefunction of 
the system. In fact, the 4— legs vertices in the curly brackets are diagonal in momentum 
space; at higher order in the loop expansion, when we contract two or more vertices of this 
type, no sum over internal loop momenta is produced, so that all higher order perturbation 
terms are suppressed by volume effects. However, we know that in the complete theory, 
the wavefunction is not factorized and loops contain all values of momentum. This suggests 
that, in order to get a finite hamiltonian, we need to introduce in the definition of our model 
some extra resummation of Feynmann diagrams, that is not automatically contained in 
this self-consistent HF approach. The only choice consistent with the cactus resummation 
performed in the two-point function by the HF scheme is the resummation of the 1-loop 
fish diagram in the four-point function. This amounts to the change from Ab to A and it is 
enough to guarantee the ultraviolet finiteness of the hamiltonian through the redefinition 



Hf) — > Hn 



A — Ab 4 

-^0 



3(A - Ab) 



At the same time the frequencies are now related to the widths {^'^_k) by 



e + M-'-?>\L'^{{vi) + {v\)), k>0 



iu'^ + 3XL-^{ipl) = mt + 3AbS 



(3.11) 



(3.12) 



Apart for 0{L~^) corrections, M plays the role of time-dependent mass for modes with 
/c 7^ 0, in the harmonic approximation. 

In this new setup the conserved energy reads 



E 



k>0 



fAb L S 



-D 



k>0 ^ ' 



(3.13) 



Since the gap-like equations ( p.l2|) are state-dependent, we have to perform the renor- 
malization first for some reference quantum state, that is for some specific collection of 
wavefunctions ^^k, as soon as and Ab are determined as functions A, ultraviolet finiteness 
will hold for the entire class of states with the same ultraviolet properties of the reference 
state. Then an obvious consistency check for our HF approximation is that this class is 
closed under time evolution. 

Rather than a single state, we choose as reference the family of gaussian states 
parametrized by the uniform expectation value {(f){x)) = L~^^'^{lpq) = (recall that we 
have {(fk) = when k ^ hj assumption) and such that the HF energy E is as small as 
possible for fixed 0. Then, apart from a translation by L^^'^cf) on yjo, these gaussian ipk are 
ground state eigenf unctions of the harmonic Hamiltonians obtained from Hk by dropping 
the quartic terms. Because of the fc^ in the frequencies we expect these gaussian states 
to dominate in the ultraviolet limit also at finite volume (as discussed above they should 
dominate in the infinite-volume limit for any k ^ 0). Moreover, since now 

{^l)=L^^' + ^, {^lk) = 7^. MO (3.14) 



2uJn 



2uk 
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the relation ( p.l2| ) between frequencies and widths turn into the single gap equation 



M-™^, + 3A,|^'^ + ^E^-,p^| (3.15) 



fixing the self-consistent value of M as a function of 0. It should be stressed that ( p.l2| ) 



turns through eq. (|3.14|) into the gap equation only because of the requirement of energy 



minimization. Generic il^k, regarded as initial conditions for the Schroedinger equations 
(p.5|), are in principle not subject to any gap equation. 

The treatment now follows closely that in the large approximation the only dif- 
ference being in the value of the coupling, now three times larger. In fact, in case of 0{N) 
symmetry, the quantum fluctuations over a given background {4>{x)) = (f) decompose for 
each k into one longitudinal mode, parallel to 0, and — 1 transverse modes orthogonal 
to it; by boson combinatorics the longitudinal mode couples to (f) with strength 3Ab/A^ 
and decouple in the N oo limit, while the transverse modes couple to (p with strength 
(A^ — l)Ab/A^ — > Ab; when = 1 only the longitudinal mode is there. 

As L ^ oo, ujI —>■ k'^ + and M is exactly the physical mass gap. Hence it must 
be A— independent. At flnite L we cannot use this request to determine and Ab, since, 
unlike M, they cannot depend on the size L. At inflnite volume we obtain 



, + 3A,[^^ + /„(M^ A)] . A) . (^^Tp^ 



d^k 1 



(3.16) 



When = this equation flxes the bare mass to be 

ml = m'^- 3Ab/D(m2, A) (3.17) 

where m = M(0 = 0) may be identifled with the equilibrium physical mass of the scalar 
particles of the inflnite-volume Fock space without symmetry breaking (see below). Now, 
the coupling constant renormalization follows from the equalities 



M^ = m^ + 3Ab[0' + Id{M\ A) - /^(m^. A)] 
= + 3A 0^ + 3A [/d(M^ A) - /^(m^ A) 



(3.18) 

finite 



A 3A , 2A , , 

- = 1-— log^ 3.19 



and reads when D = ?> 



that is the standard result of the one-loop renormalization group p2|. When D = 1, that is 
a 1 + 1— dimensional quantum fleld theory, /^(M^, A) — l£){m?,K) is already flnite and the 
dimensionfull coupling constant is not renormalized, Ab = A. 

The Landau pole in Ab prevents the actual UV limit A ^ cxo. Nonetheless, neglecting 
all inverse powers of the UV cutoff when D = 3, it is possible to rewrite the gap equation 
( CT ) as 



2 



A(M) \{m] 



+ 3 02 (3.20) 
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in terms of the one-loop running coupling constant 



A(^) = A 



1 



3A fi 



m 



It is quite clear that the HF states for which the renormalization just defined is sufficient 
are all those that are gaussian-dominated in the ultraviolet, so that we have [compare to 



eq. (|3l4]) ] 



2uJk 



er^A\ A 



oo 



(3.21) 



If this property holds at a certain time, then it should hold at all times, since the 
Schroedinger equations ( p.5| ) are indeed dominated by the quadratic term for large Uk and 
ul k'^ + const + 0{k^^) as evident from eq. ( |3.9|) . Thus this class of states is indeed closed 
under time evolution and the parameterizations ( |3.17| ) and (|3.19| ) make our tdHF approx- 
imation ultraviolet finite. Notice that the requirement (|3.21 ) effectively always imposes a 
gap equation similar to eq. ( p. 151 ) in the deep ultraviolet. 

Another simple check of the self-consistency of our approach, including the change in 
selected places from Ab to A, as discussed above, follows from the energy calculation for 
the gaussian states with {(p{x)) = (j) introduced above. Using eq. ( p. 4]) and the standard 
replacement of sums by integrals in the infinite volume limit, we find 



£:(0)= hm ^ = i02(M2-A02 



4>^ + Id{M\A) 



where M = M[(j)) depends on through the gap equation (|3.18|) . The explicit calculation 



of the integrals involved shows that the energy density difference £{4>) — £{0) [which for 
unbroken symmetry is nothing but the effective potential V^ff(0)], is indeed finite in the 
limit A — > oo, as required by a correct renormalization scheme. Notice that the finiteness of 
the energy density difference can be shown also by a simpler and more elegant argument, as 
presented below in section [illlj| . This check would fail instead when D = 3 if only the bare 
coupling constant Ab would appear in the last formula. 

The tdHF approximation derived above represents a huge simplification with respect to 
the original problem, but its exact solution still poses itself as a considerable challenge. As a 
matter of fact, a numerical approach is perfectly possible within the capabilities of modern 
computers, provided the number of equations ( |375| ) is kept in the range of few thousands. 
As will become clear later on, even this numerical workout will turn out not to be really 
necessary in the form just alluded to, at least for the purposes of this paper. 



A. On symmetry breaking 

Quite obviously, in a finite volume and with a UV cutoff there cannot be any symmetry 
breaking, since the ground state is necessarily unique and symmetric when the number of 



degrees of freedom is finite However, we may handily envisage the situation which 



would imply symmetry breaking when the volume diverges. 
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Let us first consider the case that we would call of unbroken symmetry. In this case 
the HF ground state is very close to the member with = of the family of gaussian 
states introduced before. The difference is entirely due to the quartic terms in H},. This 
correction vanish when L oo, since all wavefunctions ipk have L— independent widths, 
so that one directly obtains the symmetric vacuum state with all the right properties of 
the vacuum (translation invariance, uniqueness, etc.) upon which a standard scalar massive 
particle Fock space can be based. The HF approximation then turns out to be equivalent 
to the resummation of all "cactus diagrams" for the particle self-energy In a finite 



volume, the crucial property of this symmetric vacuum is that all frequencies uj\ are strictly 
positive. The generalization to non-equilibrium initial states with 7^ is rather trivial: 
it amounts to a shift by L^/'^cj) on '?/'o(v^o)- In the limit L — 00 we should express as 
a function of ^ = L~^^'^Lpo so that, |'?/'o(OP ~^ ^(^ ~ 0); while all other wavefunctions ipk 
will reconstruct the gaussian wavefunctional corresponding to the vacuum |0, M) of a free 
massive scalar theory whose mass M = M{(f)) solves the gap equation ( |3.18| ). The absence 
of ipQ in |0, M) is irrelevant in the infinite volume limit, since {ip^ = L^(f)'^+ terms of order 
L°. The effective potential Ves{4>) = S{(i>) — S{0), where £{4>) is the lowest energy density 
at fixed (j) and infinite volume, is manifestly a convex function with a unique minimum in 
= 0. 

Now let us consider a different situation in which one or more of the are negative. 
Quite evidently, this might happen only for k small enough, due to the /c^ in the gap equation 
[thus eq. (|3.21|) remains valid and the ultraviolet renormalization is the same as for unbroken 
symmetry]. Actually we assume here that only uJq < 0, postponing the general analysis. 
Now the quartic term in Hq cannot be neglected as L — 00, since in the ground state ipQ is 
symmetrically concentrated around the two minima of the potential IcUqV^o + iT^ V^O' that is 
ifo = ±{~uJqL^ / Xy/"^. If we scale ipo as ipo = L^^'^i then Hq becomes 

so that the larger L grows the narrower ipo{C) becomes around the two minima ^ = 
±(— cJo/A)^/^. In particular (,^^) —^l/X when L 00 and (v^q) ^ L^{^'^). More- 
over, the energy gap between the ground state of Hq and its first, odd excited state as well 
as difference between the relative probability distributions for ^ vanish exponentially fast in 
the volume L^. 

Since by hypothesis all ul with k ^ are strictly positive, the ground state ipk with 
k ^ are asymptotically gaussian when L — 00 and the relations ( |3.12|) tend to the form 



u;l = e + M^ = e + m^ 

= -2ujI = ml + 3KiL-^{ipl) + Sq) = ml + SAbC^o^ + 3KlDim\ A)] 

This implies the identification cUq = — and the bare mass parameterization 

ml=(l- |Ab/A) - SX^Ioim^ A) (3.23) 



characteristic of a negative cJq [compare to eq. ( |3.17| )], with m the physical equilibrium 



mass of the scalar particle, as in the unbroken symmetry case. The coupling constant 



12 



renormalization is the same as in eq. ( |3.19| ) as may be verified by generalizing to the 



minimum energy states with given field expectation value 0; this minimum energy is nothing 
but the HF effective potential Vg§^(0), that the effective potential in this non-gaussian HF 
approximation; of course, since ipo is no longer asymptotically gaussian, we cannot simply 
shift it by L^^'^cj) but, due to the concentration of ^/jq on classical minima as L ^ oo, one 
readily finds that Ves{4>) is the convex envelope of the classical potential, that is its Maxwell 
construction. Hence we find 




< ~Jl 

2 




and the gap equation for the 0— dependent mass M can be written, in terms of the step 
function B and the extremal ground state field expectation value v = m/v^2A, 

= m^ + 3Ab(0^ - v^) 6(02 _ ^ \Id{M\ A) - loim^ A)] (3.24) 



We see that the specific bare mass parameterization ( p.23|) guarantees the non 



renormalization of the tree-level relation v"^ = jlX ensuing from the typical symmetry 
breaking classical potential V(0) = |A(02— ■y2)2_ With the same finite part prescription as in 
eq. ( |3.18| ), the gap equation ( |3.24| ) leads to the standard coupling constant renormalization 
( Fl9| ) when D = 3. 

In terms of the probability distributions 1^^0(0 P scaled amplitude ^ = L^^^'^lpq^ 

the Maxwell construction corresponds to the limiting form 

L^oo I - 0) , 0^ > t;"^ 

On the other hand, if is indeed the only negative squared frequency, the k part of 
this minimum energy state with arbitrary = {4>{x)) is better and better approximated 
as A —s> 00 by the same gaussian state |0, M) of the unbroken symmetry state. Only the 
effective mass M has a different dependence M{(j)), as given by the gap equation ( p. 24] ) 
proper of broken symmetry. 

At infinite volume we may write 



2VF+M2 



where C(0) = 0^ in case of unbroken symmetry (that is cJq > 0), while C(0) = max(f 2, 02) 
when Uq < 0. This corresponds to the field correlation in space 

(0(x)0(y)) = / (0j(v'De^'-(^-^) = m + An{x - y, M) 

where ^^{x — y,M) is the massive free field equal-time two points function in D space 
dimensions, with self-consistent mass M. The requirement of clustering 



(0(x)0(i/)) (0(x))2 = V 



2 
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contradicts the infinite volume limit of 



k 

except at the two extremal points = ±v. In fact we know that the L oo limit of the finite 
volume states with 0^ < ^2 violate clustering, because the two peaks of -00(0 have vanishing 
overlap in the limit and the first excited state becomes degenerate with the vacuum: this 
implies that the relative Hilbert space splits into two orthogonal Fock sectors each exhibiting 
symmetry breaking, {4>{x)) = ±v, and corresponding to the two independent equal weight 
linear combinations of the two degenerate vacuum states. The true vacuum is either one 
of these symmetry broken states. Since the two Fock sectors are not only orthogonal, but 
also superselected (no local observable interpolates between them), linear combinations of 
any pair of vectors from the two sectors are not distinguishable from mixtures of states and 
clustering cannot hold in non-pure phases. It is perhaps worth noticing also that the Maxwell 
construction for the effective potential, in the infinite volume limit, is just a straightforward 
manifestation of this fact and holds true, as such, beyond the HF approximation. 

To further clarify this point and in view of subsequent applications, let us consider the 
probability distribution for the smeared field 0/ = / d^x (j){x)f{x), where 



k 



is a smooth real function with / d^x f{x) = 1 {i.e. Jq = 1) localized around the origin (which 
is good as any other point owing to translation invariance). Neglecting in the infinite volume 
limit the quartic corrections for all modes with A; 7^ 0, so that the corresponding ground 
state wavefunctions are asymptotically gaussian, this probability distribution evaluates to 

FT{u«j>f<u + du)= / rfe|^o(OPexp J ^ ^> ^ 

(ZTTZjf )-^'^ J-00 

where 



k^O 



d^k f{kf 



(27r)^ 2VA;2 + m? 



In the unbroken symmetry case we have |'0o(OP ~ — 0) as L 00, while the limiting 
form ( |3.25|) holds for broken symmetry. Thus we obtain 



PT{u<(f)f <u + du) = pf{u — (p) du , p/(u) = (27rSj) ' exp 
for unbroken symmetry and 



Pr(u <(f)f<u + du) 



^{1 + (J)/ v) p f{u — v) du + ^{1 — (f)/v) p f{u + v) du , 02 < t;2 
\pf{u — (j)) du , 02 > y2 



for broken symmetry. Notice that the momentum integration in the expression for S f needs 
no longer an ultraviolet cutoff; of course in the limit of delta-like test function /(x), 
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diverges and Pf{u) flattens down to zero. The important observation is that Pr('u < 0/ < 
u + du) has always a single peak centered in m = for unbroken symmetry, while for broken 
symmetry it shows two peaks for 0^ < v"^ and small enough. For instance, if = 0, then 
there are two peaks for Hf < v"^ [implying that f{k) has a significant support only up to 
wavevector k of order v, when D = 3, or mexp(const f^) when D = 1]. 

To end the discussion on symmetry breaking, we may now verify the validity of the 
assumption that only Uq is negative. In fact, to any squared frequency ujI (with k ^ 0) that 
stays strictly negative as L — > oo there corresponds a wavefunction ipk that concentrates on 
(fl + ip"^!^ = —ujIL^I a ; then eqs. ( |3.12| ) implies —IujI = fc^ + for such frequencies, while 
ujI = k"^ + rn^ for all frequencies with positive squares; if there is a macroscopic number of 
negative uf. (that is a number of order L^), then the expression for cUq in eq. ( p.l2|) will 
contain a positive term of order in the r.h.s., clearly incompatible with the requirements 
that cJq < and be independent of L; if the number of negative ul is not macroscopic, 
then the largest wavevector with a negative squared frequency tends to zero as L — >• cx) (the 
negative ul clearly pile in the infrared) and the situation is equivalent, if not identical, to 
that discussed above with only cJq < 0. 



B. Out— of— equilibrium dynamics 

We considered above the lowest energy states with a predefinite uniform field expectation 
value, {4>{x)) = 0, and established how they drastically simplify in the infinite volume limit. 
For generic these states are not stationary and will evolve in time. By hypothesis ipk is 
the ground state eigenfunction of when k > 0, and therefore \ipk\'^ would be stationary 
for constant Uk, but ipQ is not an eigenfunction of Hq unless = 0. As soon as \'ipk\'^ starts 
changing, {ipD changes and so do all frequencies Uk which are coupled to it by the eqs. 
( P.12| ). Thus the change propagates to all wavefunctions. The difficult task of studying 
this dynamics can be simplified with the following scheme, that we might call gaussian 
approximation. We first describe it and discuss its validity later on. 

Let us assume the usual gaussian form for the initial state [see eq. (|3.14| ) and the 
discussion following it]. We know that it is a good approximation to the lowest energy state 
with given (ipo) for unbroken symmetry, while it fails to be so for broken symmetry, only as 
far as ipo is concerned, unless 0^ > f^. At any rate this is an acceptable initial state: the 
question is about its time evolution. Suppose we adopt the harmonic approximation for all 
Hk with A; > by dropping the quartic term. This approximation will turn out to be valid 
only if the width of ipk do not grow up to the order (by symmetry the center will stay 
in the origin). In practice we are now dealing with a collection of harmonic oscillators with 
time-dependent frequencies and the treatment is quite elementary: consider the simplest 
example of one quantum degree of freedom described by the gaussian wavefunction 



1/1 .s 



2 V2(t2 a 



I- q 



where s and a are time-dependent. If the dynamics is determined by the time-dependent 
harmonic hamiltonian \[—d'l + a;(t)^g^], then the Schroedinger equation is solved exactly 
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provided that s and a satisfy the classical Hamilton equations 



2 1 



It is not difficult to trace the "centrifugal" force (4(t)~^ which prevents the vanishing of a 



to Heisenberg uncertainty principle HJTS 



The extension to our case with many degrees of freedom is straightforward and we find 
the following system of equations 

" ' ~dt^ " "^'^ ^'^ ^ 4^ ' 

coupled in a mean-field way by the relations ( |3.12| ), which now read 
ul = k^ + M^- 6A L-^al , k>0 

M^ = ml + 3A, {L-^{^1) + So) , ^o = ^T.^l ^^'^^^ 

This stage of a truly quantum zero-mode and classical modes with A; > does not appear 
fully consistent, since for large volumes some type of classical or gaussian approximation 
should be considered for ipo too. We may proceed in two (soon to be proven equivalent) 
ways: 

1. We shift Lfo = L^/'^cf) + 770 and then deal with the quantum mode ?7o in the gaussian 
approximation, taking into account that we must have (r/o) = at all times. This is 
most easily accomplished in the Heisenberg picture rather than in the Schroedinger 
one adopted above. In any case we find that the quantum dynamics of y^o is equivalent 
to the classical dynamics of and o"o = iVo)^^'^ described by the ordinary differential 
equations 



^ = -..o0-A0 , ^ = -0-0 + ^ (3.28) 







where = - 3A L-^{ifl) and {(fl) = L^(j)^ + al 



We rescale <^o = L^^"^^ right away, so that Hq takes the form of eq. (|3.22|) . Then L 



00 is the classical limit such that ipoiO concentrates on ^ = which evolves according 
to the first of the classical equations in ( p.28| ). Since now there is no width associated 



to the zero-mode, is coupled only to the widths ak with 7^ by cUq = — 3A0^, 
while M'^ = ml + 3Ab(02 + Sq). 

It is quite evident that these two approaches are completely equivalent in the infinite 
volume limit, and both are good approximation to the original tdHF Schroedinger equations, 
at least provided that ctq stays such that L'^ctq vanishes in the limit for any time. In this 
case we have the evolution equations 



i^^^2Xf--M^-)i, i^.-(,= + M>,.jL (3.29) 
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mean-field coupled by the L — oo limit of eqs. ( p.27] ), namely 

= m^ + 3Ab + S - lD{m^, A) 
for unbroken symmetry [that is as in eq. ( p.l7| )] or 



rn? + 3Ab 



12 2 
9 — V 



m 



(3.30) 



(3.31) 



for broken symmetry [that is as in eq. (|3.23|) ]. In any case we define 

„ 1 ^ n r d^k 



as the sum, or integral, over all microscopic gaussian widths [N.B.:this definition differs 
from that given before in eq. ( |3.8| ) by the classical term 0^]. Remarkably, the equations 
of motion ( |3.29|) are completely independent of the ultraviolet cut-off and this is a direct 
consequence of the substitution ( |3.11| ). Had we kept the bare coupling constant everywhere 
in the expression (|3.10|) , we would now have Ab also in front of the (f)^ in the r.h.s. of the 
first of the two equations ( p.29| ) [cfr., for instance, ref. |1[]]. 

The conserved HF energy (density) corresponding to these equations of motion reads 



S = T + V, T 



E 



V 



2 ^2 1 



+ >b( 



+ S) + |Ab(02 + S)2 - iA0^ 



(3.32) 



Up to additive constants and terms vanishing in the infinite volume limit, this expression 
agrees with the general HF energy of eq. (|3.13|) for gaussian wavefunctions. It holds both 
for unbroken and broken symmetry, the only difference being in the parameterization of the 
bare mass in terms of UV cutoff and physical mass, eqs. ( p.l7|) and ( p.23|) . The similarity 
to the energy functional of the large approach is evident; the only difference, apart from 
the obvious fact that is a single scalar rather than a 0{n) vector, is in the mean-field 
coupling o"fc-0 and o"fc-S, due to different coupling strength of transverse and longitudinal 
modes (cfr. ref. ||^). 

This difference between the HF approach for discrete symmetry {i.e N = 1) and the 
large N method for the continuous 0(A^)-symmetry is not very relevant if the symmetry is 
unbroken [it does imply however a significantly slower dissipation to the modes of the back- 
ground energy density]. On the other hand it has a drastic consequence on the equilibrium 
properties and on the out-of-equilibrium dynamics in case of broken symmetry (see below), 
since massless Goldstone bosons appear in the large N approach, while the HF treatment 
of the discrete symmetry case must exhibits a mass also in the broken symmetry phase. 

The analysis of physically viable initial conditions proceeds exactly as in the large N 
approach [^ and will not be repeated here, except for an important observation in case of 
broken symmetry. The formal energy minimization w.r.t. cr^ at fixed (p leads again to eqs. 



CTfc = 



(3.33) 
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and again these are acceptable initial conditions only if the gap equation that follows from 
eq. (|3.31|) in the L — oo limit, namely 

= m^ + 3Ab -v'^ + Id{M\ A) - Inim^, A)] (3.34) 

admits a nonnegative, physical solution for M^. Notice that there is no step function in eq. 
( p.34| ), unlike the static case of eq. ( |3.24| ), because ctq was assumed to be microscopic, so 



that the infinite volume al has no delta-like singularity in = 0. Hence M = m solves 
eq. (|3.34|) only at the extremal points = ±f , while it was the solution of the static gap 



equation ( p. 24 ) throughout the Maxwell region —v < 4> < v. The important observation is 



that eq. (|3.34|) admits a positive solution for also within the Maxwell region. In fact it 

3-power corrections 



can be written, neglecting as usual the inverse-power corrections in the UV cutoff 



A(M) A 

and there exists indeed a positive solution smoothly connected to the ground state, 
0^ = and = m^, whenever 0^ > v'^/3. The two intervals v"^ > cjP' > v"^ /?> correspond 
indeed to the metastabihty regions, while 0^ < v"^ /?> is the spinodal region, associated to a 
classical potential proportional to (0^ — v'^Y- This is another effect of the different coupling 
of transverse and longitudinal modes: in the large N approach there are no metastabihty 
regions and the spinodal region coincides with the Maxwell one. As in the large N approach 
in the spinodal interval there is no energy minimization possible, at fixed background and 
for microscopic widths, so that a modified form of the gap equation 



= + 3A 



(3.36) 



should be applied to determine ultraviolet-finite initial conditions. 

The main question now is: how will the gaussian widths at grow with time, and in 
particular how will ctq grow in case of method 1 above, when we start from initial conditions 
where all widths are microscopic? For the gaussian approximation to remain valid through 
time, all a^, and in particular ctq, must at least not become macroscopic. In fact we have 
already positively answered this question in the large N approach |§] and the HF equations 
( p.29| ) do not differ so much to expect the contrary now. In particular, if we consider the 

special initial condition = = 0, the dynamics of the widths is identical to that in the 
large N approach, apart from the rescaling by a factor of three of the coupling constant. 

In fact, if we look at the time evolution of the zero-mode amplitude ctq [see Fig. we 
can see the presence of the time-scale tl at which finite volume effects start to manifest. 
The time scale turns out to be proportional to the linear size of the box L and its 
presence prevents ctq from growing to macroscopic values. Thus our HF approximation 
confirms the large approach in the following sense: even if one considers in the variational 
ansatz the possibility of non-gaussian wavefunctionals, the time evolution from gaussian and 
microscopic initial conditions is effectively restricted for large volumes to non-macroscopic 
gaussians. 

Strictly speaking, however, this might well not be enough, since the infrared fluctuations 
do grow beyond the microscopic size to become of order L [see Fig. ^, where the evolution 
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of the mode with momentum k = 2ti/L is plotted]. Then the quartic term in the low— /c 
Hamiltonians Hj. is of order L and therefore it is not negligible by itself in the L ^ oo limit, 
but only when compared to the quadratic term, which for a fixed ul of order 1 would be 
of order L^. But we know that, when = 0, after the spinodal time and before the r^, 
the effective squared mass oscillates around zero with amplitude decreasing as and 
a frequency fixed by the largest spinodal wavevector. In practice it is "zero on average" 
and this reflect itself in the average linear growth of the zero-mode fluctuations and, more 
generally, in the average harmonic motion of the other widths with non-zero wavevectors. In 
particular the modes with small wavevectors of order feel an average harmonic potential 
with ujI of order L^^. This completely compensate the amplitude of the mode itself, so that 
the quadratic term in the low— A; Hamiltonians Hk is of order L^, much smaller than the 
quartic term that was neglected beforehand in the gaussians approximation. Clearly the 
approximation itself no longer appears fully justified and a more delicate analysis is required. 
We intend to return on this issue in a future work, restricting ourselves in the next section 
to the gaussians approximation. 



IV. LATE TIME EVOLUTION AND DYNAMICAL MAXWELL 

CONSTRUCTION 

By definition, the gaussian approximation of the effective potential Vcs{4') coincides with 
the infinite-volume limit of the potential energy V(0, {ck}) of eq. ( p.32| ) when the widths are 
of the 0— dependent, energy-minimizing form (|3.33| ) with the gap equation for admitting 
a nonnegative solution. As we have seen, this holds true in the unbroken symmetry case for 
any value of the background (f), so that the gaussian Ves is identical to the HF one, since 
all wavefunctions ipk are asymptotically gaussians as L — oo. In the presence of symmetry 
breaking instead, this agreement holds true only for 0^ > v"^; for < 0^ < w ^ the gaussian 
Ves exists but is larger than the HF potential Vg^^, which is already fiat. In fact, for any 
02 > v'^/3, we may write the gaussian Ves as 

Kff(0) = KfT(-0) = V,siv) + / duu[M{u)^ - 2Am2] 

where M{u)^ solves the gap equation (|3.35|) , namely M{u)^ = \{M{u)){3u^ — v'^). In each 
of the two disjoint regions of definition this potential is smooth and convex, with unique 
minima in +v and —v, respectively. These appear therefore as regions of metastability 
(states which are only locally stable in the presence of a suitable uniform external source). 
The HF effective potential is identical for 0^ > v'^, while it takes the constant value Veff{v) 
throughout the internal region 0^ < f^. It is based on truly stable (not only metastable) 
states. The gaussian V^fj cannot be defined in the spinodal region 0^ < v'^/3, where the gap 
equation does not admit a nonnegative solution in the physical region far away from the 
Landau pole. 

Let us first compare this HF situation with that of large There the different 

coupling of the transverse modes, three time smaller than the HF longitudinal coupling, has 



two main consequences at the static level: the gap equation similar to (|3.35|) does not admit 



nonnegative solutions for 0^ < f^, so that the spinodal region coincides with the region in 
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which the effective potential is flat, and the physical mass vanishes. The out-of-equilibrium 
counterpart of this is the dynamical Maxwell construction: when the initial conditions are 
such that (f) has a limit for t — >■ oo, the set of all possible asymptotic values exactly covers 
the flatness region (and the effective mass vanishes in the limit). In practice this means that 
101 is not the true dynamical order parameter, whose large time limit coincides with f , the 
equilibrium fleld expectation value in a pure phase. Rather, one should consider as order 
parameter the renormalized local (squared) width 

.i,„ = + = + 

Af^oo N A 

where the last equality follows from the deflnition itself of the effective mass M (see ref. 

— 2 

01). Since M vanishes as t — oo when tends to a limit within the flatness region, 
we flnd the renormalized local width tends to the correct value v which characterizes the 
broken symmetry phase, that is the bottom of the classical potential. We may say that 
the spinodal region, perturbatively unstable, at the non-perturbative level corresponds to 
metastable states, all reachable through the asymptotic time evolution with a vanishing 
effective mass. 

In the HF approximation, where at the static level the spinodal region cfp' < is 
smaller than the flatness region 0^ < f^, the situation is rather different. Our numerical 
solution shows that, (f) oscillates around a certain value 0oo with an amplitude that decreases 
very slowly. As in large A^, the asymptotic value 0oo depends on the initial value 0(0). But, 
if the background starts with zero velocity from a non-zero value inside the spinodal 
interval, then it always leaves this region and eventually oscillates around a point between 
the spinodal point v/ -\/3 and the minimum of the tree level potential v (see Fig.s ^ and 
In other words, if we start with a in the interval [— f , f], except the origin, we end up with 
a in the restricted interval [— f , — f/-\/3] U [f/\/3, f]. The spinodal region is completely 
forbidden for the late time evolution of the mean fleld, as is expected for an unstable region. 
We stress that we are dealing with true flxed points of the asymptotic evolution since the 
force term on the mean fleld [cfr. eq. ( |3.29| ), / = (2A(/)^ — M^) 0] does vanish in the limit. 
In facts its time average / = f{t)dt/T tends to zero as T grows and its mean squared 
fluctuations around / decreases towards zero, although very slowly (see Fig.s || and H). 
Moreover, for = 1 the order parameter reads as t — > oo 

(0(a:)2)K = 0^ + Sk = y + ^ , Sr = ^^ (4.1) 

where the last equality is valid for the asymptotic values and follows from the vanishing of 
the force term /. From the last formula we see that when = at the beginning, and 
then at all times, the renormalized back-reaction tends to not v"^. It "stops at the 

spinodal line". The same picture applies for a long time, all during the "slow rolling down" 
(see section to evolutions that start close enough to = 0. This fact is at the basis of 
the so-called spinodal inflation p^ |. 



In any case, the dynamical Maxwell construction, either complete or partial, poses an 
interesting question by itself. In fact it is not at all trivial that the effective potential, in any 
of the approximation previously discussed, does bear relevance on the asymptotic behavior 
of the inflnite-volume system whenever a flxed point is approached. Strictly speaking in 
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fact, even in such a special case it is not directly related to the dynamics, since it is obtained 
from a static minimization of the total energy at fixed mean field, while the energy is not 
at its minimum at the initial time and is exactly conserved in the evolution. On the other 
hand, if a solution of the equations of motion (|3.29|) exists in which the background tends 
to a constant (poo as t oo, one might expect that the effective action (which however 
is nonlocal in time) somehow reduces to a (infinite) multiple of the effective potential, so 
that should be an extremal of the effective potential. This is still an open question that 
deserves further analytic studies and numerical confirmation. 

It is worth noticing also that when the field starts very close to the top of the potential 
hill, it remains there for a very long time and evolves through a very slow rolling down, 
before beginning a damped oscillatory motion around a point in the metastability region. 
During the slow roll period, oscillates around zero with decreasing amplitude and the 
"phenomenology" is very similar to the evolution from symmetric initial conditions, as can 
be seen comparing Fig.s |^ and |^. Fig. |^ shows the evolution of the zero mode amplitude 
in case of a very slow rolling down. In such a case, after a very short (compared to the 
time scale of the figure) period of exponential growth (the spinodal time), the quantum 
fluctuations start an almost linear growth, very similar to the evolution starting from a 
completely symmetric initial state. This, obviously, corresponds to the vanishing of the 
effective mass. In the meanwhile, (f) keeps growing and rolling down the potential hill with 
increasing speed towards the minimum of the classical potential, eventually entering the 
metastable region. At that time, the effective mass starts to increase again and the zero 
mode stops its linear growth, turns down and enters a phase of "wild" evolution. This time 
scale, let us call it Tgrd, depends on the initial value of the condensate: the smaller (pit = 0) 

1 /2 

is, the longer Tgrd will be. We find numerically that Tgrd oc [(p{t = 0)j 

If we now study the dynamics in finite volume, starting from condensates different from 
zero, we will find a competition between Tsrd and tl, the time scale characteristic of the 
finite volume effects, that is proportional to the linear size of the box we put the system 
in. Fig. |l^ shows clearly that when L/2'k = 100 and (f) = 10~^, we have Tgrd ~ tl- In any 
case, either one or the other effect will prevent the zero mode amplitude from growing to 
macroscopic values for any initial condition we may start with. 

It should be noted, also, that the presence of the time scale Tsrd does not solve the 
internal inconsistency of the gaussian approximation described above in section [II B . In 
fact, for any fixed value L for the linear size of the system, we can find a whole interval of 
initial conditions for the mean field, which leave enough time to the fluctuations for growing 
to order L, much before the field itself had rolled down towards one of the minima of the 
classical potential. For those particular evolutions, we would need to consider the quartic 
terms in the hamiltonians that the gaussian approximation neglects, as already explained. 

In addition, there will be also initial conditions for which tl > Tsrd- In that case, the 
effective mass soon starts oscillating around positive values and it is reasonable to think that 
it will take a much longer time than t^ for the finite volume effects to manifest. In [^] we 
have interpreted the proportionality between tl and L as an auto interference effect (due 
to periodic boundary conditions) suffered by a Goldstone boson wave, traveling at speed of 
light, at the moment it reaches the borders of the cubic box. Here, the massless wave we 
have in the early phase of the evolution, rapidly acquires a positive mass, as soon as the 
condensate rolls down; this decelerates the wave's propagation and delays the onset of finite 



21 



volume effects. The gaussian approximation appears to be fully consistent when we limit 
ourselves to the evolution of these particular configurations. 



V. NUMERICAL ANALYSIS 

We discuss in this section the asymptotic behavior of the dynamical evolution as it turns 
out from our numerical results in the gaussian approximation. 

Let us begin with the precise form of the evolution equations for the field background 
and the quantum mode widths, as described in sections |111 B| (cfr eq. ( |3.29D . 







,2 

L h 
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(5.1) 



where the index n labels the discrete set of values used to perform the sum (finite volume) or 
the integral (infinite volume) over momenta in the quantum back-reaction S, while M^(t) is 
defined by the eq. (|3.30|) in case of unbroken symmetry and by eq. (|3.31|) in case of broken 
symmetry. The back-reaction S reads, in the notations of this appendix 



Y.9n 

?i=0 
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where Qn is the appropriate "degeneracy" factor and M is the number of modes with distinct 
dynamics. Technically it is simpler to treat an equivalent set of equations, which are formally 
linear and do not contain the singular Heisenberg term oc cr~^. This is done by introducing 
the complex mode amplitudes Zn = anexp{i9n), where the phases 9n satisfy a'^On = 1. Then 
we find a discrete version of the equations studied for instance in ref [H, namely 



subject to the Wronskian condition 







1 



Af 



(5.2) 



n=0 



One realizes that the Heisenberg term in corresponds to the centrifugal potential for the 
motion in the complex plane of Zn- Looking at the figs. ^ or ^, we can see that the motions 
of the quantum modes correspond qualitatively to orbits with very large eccentricities. In 
fact, there are istants in which (T„ is very little and the angular velocity 6'„ is very large. 
This is the technical reason for preferring the equations in the form ( |5.2| ). 

To solve these evolution equations, we have to choose suitable initial conditions respecting 
the Wronskian condition. In case of unbroken symmetry, the requirement of minimum energy 
for the quantum fluctuations leads to the massive particle spectrum: 



1 



V2n; 



dt 



(0) 




where fi„ = ^^k'^ + M'^ifi) and the initial squared effective mass M^(t 
determined self-consistently, by means of its definition (|3.3CI|) . 



0), has to be 
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In case of broken symmetry, the gap equation is a viable mean for fixing the initial 
conditions only when lies outside the spinodal region [cfr. eq ( 3.35 )]; otherwise, the gap 
equation does not admit a positive solution for the squared effective mass and we cannot 
minimize the energy of the fluctuations. Following the discussion presented in fll B| , one 
possible choice is to set o"| = ^yp^^jj < \M^\ and then solve the corresponding gap 

equation ( p.36| ). We will call this choice the "flipped" initial condition. An other acceptable 
choice would be to solve the gap equation, setting a massless spectrum for all the spinodal 
modes but the zero mode, which is started from an arbitrary, albeit microscopic, value. This 
choice will be called the "massless" initial condition. 

Before passing to discuss the influence of different initial conditions on the results, let 
us present the asymptotic behavior we find when we choose the flipped initial condition. In 
Fig. |l^ we have plotted the asymptotic values of the mean field versus the initial values, 
for A = 0.1. All dimensionful quantities are expressed in terms of the suitable power of the 
equilibrium mass m. For example, the vev of the field is equal to a/5 in these units. First of 
all, consider the initial values for the condensate far enough from the top of the potential hill, 
say between 0(t = 0) = 0.88 and 0(t = 0) = 2.64. In that region the crosses seem to follow a 
smooth curve, that has its maximum exactly at 0oo = (the point of stable equilibrium). 
When we start from an initial condition smaller than 0(t = 0) = 0.88, the asymptotic value 
0OO is not guaranteed to be positive anymore. On the contrary, it is possible to choose 
the initial condition in such a way that the condensate will oscillate between positive and 
negative values for a while, before settling around an asymptotic value near either one or 
the other minimum, as fig ^ clearly shows. Fig.s |15| and |1^ helps to understand this 
behavior by consideration on the energy balance. Both the evolutions are such that the 
classical energy, defined as A(0^ — f^)/4, is not a monotonically decreasing function of time. 
Indeed, energy is continuously exchanged between the classical degree of freedom and the 
quantum fluctuations bath, in both directions. However, the two rates of energy exchange 
are not exactly the same and an effective dissipation of classical energy on average can be 
seen, at long time at least. Of course, this is not the case for the initial transient part 
of the evolution starting from the initial condition 0(t = 0) = 0.08; there, the condensate 
absorbs energy (on average) from the quantum fluctuations, being able to go beyond the 
top of the potential hill, towards the negative minimum. This happens because in case of 
broken symmetry, the minimization of the fluctuation energy, within microscopic gaussian 
states, is not possible for initial conditions in the spinodal region [cfr. the discussion about 
the gap equation ( 3.35|) in section IIIB ]. After a number of oscillations, the energy starts 
to flow from the condensate to the quantum bath again (on the average), constraining the 
condensate to oscillate around a value close to one of the two minima. If we look at fig. ^ 
again, we can find positive asymptotic values as well as negative ones, without a definite 
pattern, in the whole interval [0.01,0.8]. If we start with < 0(t = 0) < 0.01 we have the 
slow rolling down, already described in section |^ and the mean field oscillates around a 
positive value from the beginning, never reaching negative values. A further note is worth 
being added here. During the phase of slow rolling down, the evolution is very similar to 
a symmetric evolution starting from (j){t = 0) = 0; in that case, the dissipation mechanism 
works through the emission of (quasi-) massless particles and it is very efficient because it has 
not any perturbative threshold. If the field stays in this slow rolling down phase for a time 
long enough, it will not be able to absorb the sufficient energy to pass to the other side ever 
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again and it will be confined in the positive valley for ever. Evidently, when 0(t = 0) > 0.01 
this dissipative process might not be so efficient to prevent the mean field from sampling 
also the other valley. 

Which one of the two valleys will be chosen by the condensate is a matter of initial 
conditions and it is very dependent from the energy stored in the initial state, as is shown in 
fig. 0, where two evolutions are compared, starting from the same value for the condensate, 
but with the two initial conditions, "flipped" and massless, for the quantum fluctuations. 



VI. CONCLUSIONS AND PERSPECTIVES 

In this work we have extended the standard time dependent Hartree-Fock approximation 
p!3| for the 0'^ QFT, to include some non-gaussian features of the complete theory. We have 



presented a rather detailed study of the dynamical evolution out of equilibrium, in finite 
volume (a cubic box of size L in 3D), as well as in infinite volume. For comparison, we 
have also analyzed some static characteristics of the theory both in unbroken and broken 
symmetry phases. 

By means of a proper substitution of the bare coupling constant with the renormalized 
coupling constant (fully justified by diagrammatic consideration), we have been able to 
obtain equations of motion completely independent of the ultraviolet cut-off (apart from a 
slight dependence on inverse powers, that is, however, ineluctable because of the Landau 
pole). We have described in detail the shape of the ground state, showing how a broken 
symmetry scenario can be recovered from the quantum mechanical model, when the volume 
diverges. 

Moreover, we have shown that, within this slightly enlarged tdHF approach that allows 
for non-gaussian wavef unctions, one might recover the usual gaussian HF approximation 
in a more controlled way. In fact, studying the late time dynamics, we have confirmed 
the presence of a time scale r^, proportional to the linear size L of the box, at which 
the evolution ceases to be similar to the infinite volume one. At the same time, the low- 
lying modes amplitudes have grown to order L. The same phenomenon has been observed 
in the 0{N) model ^j. Looking at this result in the framework of our extended tdHF 
approximation, one realizes that the growth of long- wavelength fluctuations to order L in 
fact undermines the self-consistency of the gaussian HF itself. In fact, in our tdHF approach 
the initial gaussian wavefunctions are allowed to evolve into non-gaussian forms, but they 
simply do not do it in a macroscopic way, within a further harmonic approximation for 
the evolution, so that in the infinite-volume limit they are indistinguishable from gaussians 
at all times. But when is on average not or order but much less, as it happens 
for suitable initial conditions, infrared modes of order L will be dominated by the quartic 
term in our Schroedinger equations (|3.5|) , showing a possible internal inconsistency of the 
gaussians approximation. 

An other manifestation of the weakness of the HF scheme is the curious "stopping at the 
spinodal line" of the width of the gaussian quantum fluctuations, when the initial configura- 
tion does not break the symmetry. This does not happen in the large N approach because 
of different coupling of transverse mode (the only ones that survive in the N oo limit) 
with respect to the longitudinal modes of the = 1 case in the HF approach. 
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We have also described the non-trivial phenomenology of the infinite-volume late-time 
evolution in the gaussian approximation, showing how the dynamical Maxwell construction 
differs from the N = oo case. In fact, we have observed the presence of an unstable interval, 
contained in the static flat region which is forbidden as attractor of the asymptotic evolution. 
This region corresponds, more or less, to the spinodal region of the classical potential, with 
the obvious exception of the origin. In particular, we have found that the energy flux between 
the classical degree of freedom and the bath of quantum fluctuations is quite complex and 
not monotonous. In other words, since we start from initial conditions where the fluctuation 
energy is not minimal, there are special situations where enough energy is transferred from 
the bath to the condensate, pushing it beyond the top of the potential hill. 

Clearly further study, both analytical and numerical, is needed in our tdHF approach to 
better understand the dynamical evolution of quantum fluctuations in the broken symmetry 
phase coupled to the condensate. An interesting direction is the investigation of the case 
of flnite N, in order to interpolate smoothly the results for = 1 to those of the 
approach. It should be noted, in fact, that the theory with a single scalar fleld contains 
only the longitudinal mode (by deflnition), while only the transverse modes are relevant in 
the large N limit. Hence a better understanding of the coupling between longitudinal and 
transverse modes is necessary. 

In this direction, another relevant point is whether the Goldstone theorem is respected in 
the HF approximation pO|- It would be interesting also to study the dynamical realization 



of the Goldstone paradigm, namely the asymptotic vanishing of the effective mass in the 
broken symmetry phases, in different models; this issue needs further study in the 2D case 
101, where it is known that the Goldstone theorem is not valid. 
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FIGURES 





FIG. 2. Next-to-zero mode (k = 271 /L) amplitude evolution for different values of the size 
= 20, 40, 60, 80, 100, for A = 0.1 and broken symmetry, with (p = 0. 
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FIG. 3. Evolution of the background for two different initial conditions within the spinodal 
interval, in the tdHF approximation, for A = 1; (j){t = 0) = 0.1 (dotted line) and (j){t = 0) = 0.4 
(solid line). 
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FIG. 4. Evolution of for the two initial conditions of fig. |^. 
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FIG. 5. The average force f, defined as f = J f{t)dt/T, plotted vs. T, for A = 0.1 and 
(j) = 10^^ (solid line), <f> = 10~^ (dashed line) and cp = 10~^ (dotted- dashed line). 
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FIG. 6. The mean squared fluctuations of the force f, defined as J^{f{t) — f)'^dt/T , plotted vs. 
T, for the three initial conditions of fig. |^. 
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FIG. 7. The evolution of the mean value (solid line), the quantum back-reaction S (dashed 
line) and the squared effective mass (dotted- dashed line), for = at t = 0. 




FIG. 8. The evolution of the mean value (solid line), the quantum back-reaction S (dashed 
line) and the squared effective mass (dotted- dashed line), for (j) = 10"^ at t = 0, and A = 0.1. 
The field rolls down very slowly at the beginning. 
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FIG. 9. Evolution of the amplitude of the zero mode for A = 0.1 and (f) = 10 




FIG. 10. Comparison between the evolutions of the zero mode amplitude in the following two 
situations: the dashed line corresponds to a finite volume simulation with L/2tt = 100 and (/> = 0, 
while the solid line refers to the infinite volume evolution, with 4> = 10~^. Both correspond to 
A = 0.1. 
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FIG. 11. Evolution of the mean value (j) for A = 0.1 and for two different initial conditions: 
4> = 0.08 (solid line) and ^ = 0.16 (dashed line), with the "flipped" choice for the spinodal modes. 




FIG. 12. Evolution of the mean value (j) for A = 0.1, with (p{t = 0) = 0.08, and two different 
initial conditions for the quantum spinodal modes, "flipped" (solid line) and massless (dashed line). 
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FIG. 13. Asymptotic values of the mean field (p, plotted vs. initial values (j){t = 0), for A = 0.1. 
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FIG. 14. Comparison between the classical energies for the two initial conditions of Fig [7J. 



32 



0.4 - li 




FIG. 15. Evolution of the condensate (j) (solid line) and of the corresponding classical energy 
(dashed line), for (j){t = 0) = 0.16 and A = 0.1 (cfr. Fig.s and IJ). 




FIG. 16. Evolution of the condensate <j) (solid line) and of the corresponding classical energy 



(dashed line), for (t){t = 0) = 0.08 and A = 0.1 (cfr. Fig.s Wa and 11). 
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